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Abstract 

I discuss the impact of gauge-gravity duality on our understanding of two classes of systems: 
conformal quantum matter and compressible quantum matter. 

The first conformal class includes systems, such as the boson Hubbard model in two spatial 
dimensions, which display quantum critical points described by conformal field theories. Ques- 
tions associated with non-zero temperature dynamics and transport are difficult to answer using 
conventional field theoretic methods. I argue that many of these can be addressed systematically 
using gauge-gravity duality, and discuss the prospects for reliable computation of low frequency 
correlations. 

Compressible quantum matter is characterized by the smooth dependence of the charge density, 
associated with a global U(l) symmetry, upon a chemical potential. Familiar examples are solids, 
superfluids, and Fermi liquids, but there are more exotic possibilities involving deconfined phases 
of gauge fields in the presence of Fermi surfaces. I survey the compressible systems studied using 
gauge-gravity duality, and discuss their relationship to the condensed matter classification of such 
states. The gravity methods offer hope of a deeper understanding of exotic and strongly-coupled 
compressible quantum states. 



I. INTRODUCTION 



One of the remarkable developments to emerge from research in string theory in the past 
decade is the idea of gauge-gravity duality pp. This is an equivalence between a quantum 
field theory in D spacetime dimensions, and a quantum theory of gravity in D + 1 spacetime 
dimensions. The D-dimensional theory does not have a gravitational force, and is to be 
viewed as a 'hologram' of the D + 1 dimensional theory. This is a true equivalence between 
theories, and not a projection to a restricted portion of Hilbert space: the number of degrees 
of freedom of the D-dimensional theory equal those of the (D + l)-dimensional theory with 
gravity. It is thus in the spirit of early ideas of Bekenstein and others [2H5] that in theories 
of quantum gravity, the entropy of a region of spacetime is proportional to its spatial surface 
area (and not its spatial volume). 

A powerful feature of the duality is that it is also a strong-weak coupling duality: it maps 
a strongly-coupled theory in D dimensions to a weakly-coupled theory in D + 1 dimensions, 
and vice versa. It has therefore, justifiably, raised hopes for understanding new classes of 
strongly coupled quantum field theories in D dimensions. This hope extends to researchers, 
like condensed matter physicists, who had no previous interest in theories of quantum gravity. 

The initial examples [1] of gauge-gravity duality appeared in highly supersymmetric gauge 
theories e.g. the SU(iV) Yang-Mills gauge theory with maximal supersymmetry (A/" = 4) 
in D = 4. This is dual to a string theory on a 5-dimensional space with constant negative 
curvature: anti-de Sitter space, abbreviated AdS 5 . In the low energy limit, the string theory 
reduces to a gravity theory on AdSs also with maximal supersymmetry. Clearly, neither of 
these models are of any specific interest to condensed matter physics. However, since then it 
has become clear that this is but the simplest of a general class of duality between conformally 
invariant field theories (CFTs) in D dimensions and theories of gravity in AdSa+i. Strong 
evidence for such a duality has also appeared for CFTs without any supersymmetry [61 [7] . 
As we will review in Section II A , CFTs appear at quantum phase transitions of generic 
models relevant for experiments in condensed matter and ultracold atomic gases, and some 
of these are directly in the classes for which gauge-gravity duality has been studied. 

It is perhaps useful here to make an analogy with earlier developments in the solution of 
strongly interacting many body systems. Starting with the work of Bethe in 1931 [8J, a wide 
class of solvable quantum many body systems were discovered in one spatial dimension (D = 
2). These models are 'integrable' in that they have an infinite number of conservation laws, 
and explicit wavefunctions can be written down for all eigenstates using the Bethe ansatz. 
However, integrability invariably required fine-tuning of the structure of the Hamiltonian, 
and no one expected that integrable systems would be directly applicable to any experiments. 
Nevertheless, the structure of the integrable models and their excitations taught us a great 
deal about quantum many body physics in D = 2, and these insights were crucial in the 
development of the far more general effective field theory of the 'Tomonaga-Luttinger' liquid 
[9]. The latter theory applies to generic interacting systems in D = 2, and has found 
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numerous experimental applications. 

The hope is that a similar success will eventually be achieved using gauge-gravity duality, 
for strongly interacting quantum systems in D > 2. We have the analog of the Bethe ansatz: 
a rapidly-increasing class of special models with known gravity duals. Much research [TUL 
IT8] is now directed towards generalizing the insights obtained from these solvable models 
to a general effective field theory approach to strongly interacting quantum systems using 
holography. 

I will review applications of gauge-gravity duality to two broad classes of condensed 
matter systems, which I denote "conformal quantum matter" and "compressible quantum 
matter" . 

In the first "conformal" class I include models with quantum critical points described 
by CFTs. Here the relevance of the holographic approach is evident, as the equivalence 
to dual gravity models has been explicitly demonstrated. Much has been learnt about 
such CFTs using traditional field-theoretic methods (such as the e expansion of Wilson and 
Fisher [IH])- However, there are key questions about their long-time correlations at non- 
zero temperatures (T) which cannot be described in a controlled manner by these methods. 
Numerical studies also fail for such questions because of difficulties associated with the "sign" 
problem: quantum simulations require summing over intermediate states whose weights 
are not positive real numbers. Remarkably, holography does yield useful results for such 
long-time correlations, including the specific numerical values of transport co-efficients and 
damping rates: these results will be discussed in Section [IT] 

The second class of "compressible quantum matter" will be defined more precisely in 
Section |III| But as the name indicates, these are systems which have a non-zero and finite 
compressibility at T = 0. Familiar examples of compressible states are Fermi liquids, solids, 
and superfluids. However, much modern condensed matter research has focused on the 
search for new types of compressible states, motivated primarily by the "strange metal" 
behavior of numerous correlated electron systems [20] ■ Current studies using holography 
have produced many examples of compressible quantum matter, although their place in the 
traditional condensed matter classification has not yet been fully established. These issues 
are the focus of active research, as I will describe in Section |III| 

I close this section by noting that the above classification implicitly assumes D > 2. The 
case D = 2 is special: Tomonaga-Luttinger liquids are both conformal and compressible, 
and this is partly responsible for their simplicity. I will implicitly assume D > 2 in the 
remainder of the paper, where the two classes are distinct and far more complicated, and 
have properties quite different from Tomonaga-Luttinger liquids. 
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II. CONFORMAL QUANTUM MATTER 



A. Condensed matter analysis 

Let us begin by describing the simplest model which realizes a CFT, realizable in a 
laboratory [211 E2]- Consider spinless bosons hopping on a lattice of sites, i, with short- 
range repulsive interactions. With a boson annihilation operator bi, such physics is usefully 
captured by the Hubbard Hamiltonian 



H b = -wJ2 ( b l b j + h ) h i) +^^2 n i( n i ~ !) 

(ij) * 



(2.1) 



where bi is the canonical boson annihilation operator, rii = b\bi is the boson number oper- 
ator, w is the hopping matrix element between nearest-neighbor sites, and U is the on-site 
repulsive energy between a pair of bosons. Let us assume that the average boson density is 
exactly 1 per site, and examine the ground state as a function of the dimensionless parameter 
U/w. For large U/w, the boson repulsion dominates and so the bosons avoid each other by 
localizing on separate sites. This leads to the insulating state shown in Fig. [TJ any motion 
of bosons requires placing at least two on a site, and this is suppressed by the repulsive 
energy. In the opposite limit of small U/w, we can treat the system as a nearly free Bose 
gas, and this undergoes Bose-Einstein condensation to a the superfluid state also illustrated 
in Fig. [TJ In this state, a snapshot of the wavefunction shows large number fluctuations on 
each lattice site, and so a supercurrent is able to flow without dissipation. 



L-% — 



FIG. 1: Superfluid (left) and insulating states of bosons with repulsive interactions on a square 
lattice at / = 1. 



The distinct ground states in the limit of small and large U/w implies they cannot be 
smoothly connected: they are separated by a quantum phase transition at an intermediate 
critical point. More generally we denote parameters like U/w by a generic coupling g, and 
the quantum critical point appears at g — g c . The nature of the ground state near and 
a t g — g c is weu understood, and a pedagogical description appears in the Supplementary 
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Material. For our purposes, the most important property is that the ground state and its 
low energy excitations are efficiently described by a universal quantum field theory [23J. 
By 'universal' we mean that the same field theory applies for a large class of models of 



the superfluid-insulator transition, and not just the Hubbard model in Eq. (2.1). The field 
theory is expressed in terms of a complex field ip(r,t), which is just the continuum limit of 
the boson annihilation operator bi, with r the D — 1 dimensional spatial co-ordinate, and t 
the time. The action for the field theory is (see Supplementary Material) 



S b = J d D - l rdt (-\d t ^\ 2 + v 2 \VM 2 + (g- g c M 2 + u\^) ; 



(2.2) 



note that its couplings change as a function of the coupling g. This low energy quantum 
field theory already has an 'emergent' symmetry not shared by the Hubbard model: it is 
invariant under Lorentz transformations, with the velocity v playing the role of the velocity 
of light. Here v is a sound velocity, and its value is determined by w, U, and the lattice 
spacing. 

The symmetry of the theory becomes much larger at the quantum critical point at g — g c . 
Here, as may be familiar to some readers from the theory of second order phase transitions, 
the theory is scale invariant. Specifically, the structure of the quantum correlations remain 
invariant under the rescaling transformation 

t->t/b , r^r/b (2.3) 

where b is the rescaling factor. Formally, this scale invariance arises because the quantum 
critical point is realized as a fixed point of the renormalization group transformation. Actu- 
ally, we are not done with the set of emergent symmetries. Quantum field theories which are 
Lorentz invariant and which obey certain 'hyperscaling' properties, are also invariant un- 
der conformal transformations of spacetime. These are transformations which preserve the 
Lorentzian metric of spacetime upto an overall factor, and so preserve all angles. Specifically, 
the Lorentzian metric rescales under conformal transformations 

(-v 2 dt 2 + dr 2 ) {-v 2 dt 2 + dr 2 ) /b 2 (r, t) (2.4) 

where the rescaling factor b can be a certain set of smooth functions of spacetime; the 
rescaling transformation is a special case of a conformal transformation. The critical point 
of Sb, describing the superfluid-insulator quantum phase transition, is invariant under such 
conformal transformations, and is so a CFT. We will restrict further attention to the case of 
2 spatial dimensions, with D — 3, when this CFT is strongly coupled. The CFTs for D > 4 
are essentially free field theories, and we don't need sophisticated methods to analyze them. 

As we noted earlier, many properties of the strongly coupled CFT at T = are well 
understood. The main tool is the renormalization group, and its realization in the context 
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of various analytic and numerical expansions — these are well-established methods which I 
will not dwell on here. However, all of these methods fail for certain key questions on the 
dynamics at T > 0. 

To describe these questions, we need the phase diagram of the model at T > [21], shown 
in Fig. [2} There is much interesting physics associated with the many distinct features of this 
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FIG. 2: Phase diagram of the superfiuid-msulator transition in two spatial dimensions {D = 3). 
The quantum critical point is at g = g c , T = 0. The dashed lines are crossovers, while the full line 
is a phase transition at the Kosterlitz-Thouless temperature Tkt > 0. 



phase diagram, but for now the reader is asked to focus on the distinction between the blue- 
and pink-shaded regions. In the blue-shaded regions, the physics can be described in terms 
of the familiar excitations of either the insulator or the superfluid, which are illustrated in 
Fig. [3j For the insulator, the excitations are particle or hole excitations above the background 
of the insulator with one particle per site. In contrast, for the superfluid, the excitations are 
point-like vortices in the background of the Bose condensate. A semiclassical theory of gases 
of such excitations provides an essentially complete description of the long-time correlations 
in the blue-shaded regions. 

Let us now turn to the pink-shaded region of Fig. [2j separated from the blue-shaded 
regions by crossovers indicated by the dashed lines; the reader is referred to another review 
article [20] for a detailed discussion of the location and shape of these crossover lines. The 
defining characteristic of the pink-shaded 'quantum critical' region is that its dynamics is 
controlled by the CFT and its excitations. These excitations do not have a particle-like 
interpretation, they are not amenable to an effective classical description, and they interact 
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FIG. 3: Excitations of the superfluid and the insulator. The excitations of the superfluid are 
vortices and anti-vortices, which can be viewed as point-like particles centered at the X. A single 
particle and single hole excitations of the insulator are shown on the right. 

strongly (and universally) with each other. We are now faced with the challenge of describing 
the long-time dynamics of this strongly interacting quantum critical regime. 

The traditional renormalization group methods have been applied to the problem of 
quantum critical dynamics [25], [26]. While a great deal of insight, and some analytic results 
have been obtained, the results are not expected to be qualitatively accurate. We will 
shortly turn to a description of the new holographic methods, and describe their promise in 
eventually solving this challenging problem. 

It is helpful to focus on a specific observable characterizing the quantum critical dynamics 
in the pink-shaded region of Fig. [2j Let us choose the frequency-dependent conductivity 
a{oj): we endow the bosons with a (possibly fictitious) charge q, and examine the current 
response in the presence of an 'electric' field coupling to the q charge, which is oscillating at 
a frequency lo. General arguments based on the properties of the CFT imply that [26] 




where h is Planck's constant/ (2tt), ks is Boltzmann's constant, and £ is an unknown, 
dimensionless universal function of the dimensionless ratio of frequency to temperature. This 
structure follows from the fact that in D = 3 the conductivity is a dimensionless number 
when expressed in units of Q 2 /h; and the ability of the strong and universal interactions 
between the excitations of the CFT to relax the electrical current. Each CFT in D = 3 
is characterized by its own function E(cJ). We are now faced with the basic challenge: 
accurately determine the function E(cJ) for the quantum critical point of the theory Sb in 





(2.5) 
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Eq. (2.2) describing the superfluid-insulator transition. As we will see below, even simple 
questions on the overall shape of this function remain unresolved. 

The traditional field-theoretic methods (which are based on the renormalization group), 
as well as numerical studies, do work accurately in determining the high frequency limit 
ftw ^> ksT. Here, the correlations are essentially unchanged from the ground state values, 
and so are amenable to a controlled analysis. The value of the number Eqo = £(a> — > oo) has 
been so determined [2TH21?] . and there are no obstacles (in principle) in refining this work 
to even more precise values. 

How about determining the opposing limit [30], the value of the number So = S(tu — > 0)? 
This is a number characterizing the long-time response at a non-zero temperature. So we 
have to extrapolate from our knowledge of the short-time dynamics to the long-time limit. 
For systems with excitations which are weakly interacting particles, such extrapolations 
have traditionally been carried out by the venerable Boltzmann equation, and its many 
descendants. Such methods are evidently not directly applicable to the non-particle-like 
excitations of a CFT, which also have strong interactions. 

Nevertheless, let us forge ahead and see what the Boltzmann approach teaches us about 
the qualitative shape of the function S. Examining Fig. [3j we are immediately faced with a 
choice: do we approach the critical point starting from the insulator or the superfluid? As 
we will now argue, a qualitatively different answer obtains from the two approaches. 

Let us begin by extrapolating the excitations of the insulator to the critical point. The 
insulator has particle and hole excitations, and there is no difficulty in writing down a 
Boltzmann equation for these excitations [26]. In the simplest picture, thermally excited 
quasiparticles undergo Brownian motion with mutual collisions, while drifting in the applied 
'electric' field E. The average velocity of these particles, v, therefore obeys an equation like 

dv v „ , . 

It + 7 C = qE (2 - 6) 

where r c is the mean time between collisions of the particles. If we extrapolate this picture to 
the critical point (without much justification), we expect that this time would be determined 
by the only energy scale characterizing the CFT, which is fc#T, and hence r c ~ h/{ksT). 



Solving Eq. (2.6, we then predict a "Drude" form for the frequency-dependent conductivity 



at low frequencies, with 

<«>) = v^—- ^ 

1 — lUJT r 



Combining the real part of Eq. (2.7) with our earlier considerations in the large u limit, 
we can surmise a frequency dependent conductivity of the CFT with the qualitative form 
shown in Fig. [1] [26]. We reiterate that this form implicitly builds in the structure of the 
particle-like excitations of the insulator. 

However, in principle, it should be equally valid to approach the critical point from the 
superfluid side, using its vortex-like excitations as the basic degree of freedom. In two spatial 
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FIG. 4: Expected form for the frequency-dependent conductivity of the CFT of the superfluid- 
insulator transition, using the Boltzmann picture applied to the particle excitations of the insulator 

dimensions, vortices can be specified by the location of their point-like center, and so can 
also be viewed as excitations which are 'particles'. We can, therefore, similarly write down 
a Boltzmann equation for the thermal and quantum dynamics of these vortex excitations. 
We will refrain from doing so here, but note a basic property we will need: the physical 
conductivity of the underlying boson model is equal to the resistivity of the vortex-like 
particles entering this Boltzmann treatment [31 J . This can be deduced from the fact that 
a flow of vortices in a superfluid induces a voltage in the transverse direction. From this 
we expect that the conductivity of the quantum critical point is roughly the inverse of the 
previous prediction in Fig. |4j this result for the vortex picture of transport is shown in Fig. [5[ 

We now have two qualitatively distinct predictions for the frequency-dependent function 
E, shown in Fig. [4] and [5} It is a fair statement that for the CFT of the superfluid-insulator 
insulator transition described by the boson Hubbard model, we do not know today which of 
these results, if either, is the correct one. The conventional e and vector large N expansions 
[26[ 132] both have a built-in bias towards the excitations of the insulator, and so predict a 
form qualitatively similar to that in Fig. |4j However, we are unable to say if these results 
should be believed for the physical case of interest. It is remarkable that this basic property 
of one of the simplest quantum critical point in two spatial dimensions remains unresolved. 
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FIG. 5: Expected form for the frequency-dependent conductivity of the CFT of the superfluid- 
insulator transition, using the Boltzmann picture applied to the vortex excitations of the superfluid. 
The "d.c. conductivity" of the vortex- like excitations is a$ v . 

B. Holographic analysis 

Let us now change gears, and discuss what gauge-gravity duality can teach us about the 
questions at the end of the last section. 

I will only be able to present a very short summary of the principles of gauge-gravity 
duality here; the reader is referred to the many other reviews in the literature [T0l - [T8] . 

The safest route is to begin with the solvable model, and then use physical arguments to 
generalize to a wider class of CFTs, including perhaps those of physical interest to us. The 
canonical solvable model in D = 3 spacetime dimensions is a highly supersymmetric CFT 
now known as the ABJM model [33]. In a particular 'large N y limit, the low energy limit 
of this CFT is known to map onto a supersymmetric gravity theory in D = 4 spacetime 
dimensions. We shall be interested here in the nature of correlations of currents of globally 
conserved charges (analogous to the electrical current of the boson model above) in this 
theory. Choosing a particular conserved current of the ABJM theory (it doesn't matter which 
one), its current correlations are described by an especially simple and familiar theory of 
gravity and 'electromagnetism' in D = 4: the Einstein-Maxwell (EM) theory with a negative 
cosmo logical constant, with action [15] 



Here x = (t, u, r) is a 4-dimensional spacetime co-ordinate upon which the gravity theory 
is based; r labels the 2-dimensional space upon which our CFT lives, t is time, and u 




(2.8) 
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is the new 'emergent' co-ordinate of spacetime. The coupling k is related to Newton's 
gravitational constant by k 2 = 8ttG. The gravity theory is expressed in terms of a metric 
g and its Riemann curvature scalar R, while the U(l) gauge theory has 'electromagnetic' 
flux F a t,] we use small Latin letters for the spacetime co-ordinates of the D + 1-dimensional 
spacetime. Here, and henceforth, we have set the velocity of 'light' to unity, v — 1. 

The physical significance of the gravity and the U(l) gauge field can be described after 



considering the vacuum solution of the classical equations of motion associated with Eq. 2.8 



This solution has F^ = 0, and the metric associated with 

ds i= f£y du * + (^y(_ dt 2 +dr 2j (2J) 

This is the metric of the space of uniform negative curvature, known as AdS4, and L is the 
radius of curvature. It is evidently invariant under Lorentz transformations, and also under 



scale transformations in (2.3) provided we choose 

u^ub; (2.10) 
furthermore, it is also invariant under a suitable extension of the conformal transformation in 



(2.4). This invariance is the crucial connection relating AdS4 to conformal transformations, 



and to CFTs: the group of isometries of AdS^+i is the same as the group of conformal 



transformations in D spacetime dimensions. The transformation in (2.10) also suggests a 
physical interpretation for the new co-ordinate u: it transforms like an energy /momentum 
scale, and so can be viewed as the running energy /momentum scale for the renormalization 
group [IS]. Thus the physics as u — > oo is the ultraviolet (UV) or short-distance/time 
physics, while the physics as u — > is the infrared (IR) or long-distance/time physics. The 
gravity theory on AdSo+i maintains the complete 'history' of the renormalization group 
flow in the structure of the metric as a function of u. 



We also note that the AdS space in Eq. (2.9) has a boundary as u — > oo. This boundary 



is just the D-dimensional Lorentzian spacetime upon which our CFT of interest 'lives'. 



We are now in a position to describe the role of the U(l) gauge field F^ in Eq. (2.8). This 
is the field which is 'dual' to the conserved current of the CFT. Let us label the conserved 
current J M ; we use small Greek letters to represent the components of the D-dimensional 
Lorentzian spacetime. The conductivity is related to the two-point correlator of J^. To 
enable computation of this two-point correlator, let us include a source term coupling to the 
conserved current of the CFT: 



Soft -> 5 cft - J dP^rdtK^ (2.11) 
Then a fundamental property of the gauge-gravity duality is that this source, K„, is the 



11 



limiting boundary value of the vector potential Am associated with the U(l) flux F a f, [T5] : 



Ap(r,t,u -» oo) 



(2.12) 



The complete prescription for computing the conductivity of the CFT using the holographic 
approach is as follows. Solve the equations of motion of the gravitation theory subject to 



the constraint in Eq. (2.12), for an arbitrary spacetime dependent K^(r,t). From these 
solutions, compute the functional dependence of the gravitational action on K^. This is 
precisely the functional dependence of the CFT action on K^, and so correlators of J M are 
easily obtained by taking functional derivatives of the action with respect to K^. 

Before we describe the results of such a computation, we need to discuss 2 additional 
points. 

First, the discussion above has implicitly assumed that we were at T = 0. However, 
our most difficult questions about the CFT were at T > 0, so it is essential we extend to 
non-zero temperatures. The key to this extension is to examine black hole solutions of the 



equations of motion of Eq. (2.8). These are analogous to the Schwarzschild solution, and 
take the form p21 [15] 



ds' 



where 



du 2 



/(«) 



+ 



(-f(u)dt 2 + dr 2 ) 



1 



(2.13) 



(2.14) 



in D = 3. Here R is a parameter which labels the location of the black hole horizon. 
Also, strictly speaking, this is not a black hole, but a black brane: the horizon is spatially 
infinite, and extends across the flat two-dimensional r space. As shown by Hawking [3], the 
black brane horizon must have a temperature, T, and in the present duality we can identify 
this temperature with the temperature, T, of the CFT. This Hawking temperature is most 



simply computed by analytically continuing the metric in Eq. (2.13) to imaginary time, and 



demanding that the resulting space be periodic in the imaginary time direction with period 
H/(kBT); such a computation yields 



k B T 



3hvR 
AnL 2 



(2.15) 



where we have momentarily reinserted the factor of v; this can be viewed as an equation 
which fixes the value of R. This black-brane is a powerful feature of the present holographic 
approach. It enables us to see how dissipative and relaxational processes of the CFT at non- 
zero temperature have a natural interpretation in the gravitational theory. As illustrated 
in Fig. M waves propagating in the D + 1 dimensional space get damped because they lose 



energy across the black-brane horizon: it is this damping which Eq. (2.12) eventually relates 
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to the dissipative transport co-efficients of the CFT 
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FIG. 6: Schematic of the geometry of gauge-gravity duality for conformal quantum matter. The 
gravity theory is defined on the bulk 4-dimensional spacetime (the time co-ordinate, t, is not 
shown). The CFT resides on the boundary at u — > oo. Waves falling across the black brane at 
u = R capture the dissipation of the CFT at non-zero temperature. 



Second, as we noted earlier, the action in Eq. (2.8) was obtained for the special CFT 



described by ABJM. It also describes a very wide class of related supersymmetric CFTs, 



but this does not include the CFT associated with the boson Hubbard model in Eq. (2.2) 



To extend to this wider class, we will use the spirit of effective field theory, applied to 
the holographic method as proposed recently in Ref. ESJ We can also view Eq. (2.8) as 



the simplest effective action for the metric and the U(l) gauge field, invariant under the 
underlying symmetries, in which all fields are expanded to include 2 gradients of spacetime 
co-ordinates. Let us assume this is a reasonable starting point, and extend this action to 
include 4 spacetime gradients. As we are only interested in linear response to the source 



Kp in Eq. (2.11), we can restrict this extension to include only 2 powers of Am- With these 
conditions, it turns out only one additional term is allowed, modulo allowed co-ordinate 
transformations, and the extended action is 



S = S E m + / d A x 



(2.16) 



where is the Weyl curvature tensor. Eq. ( 2. 16 ) is our final theory for a very wide class of 



CFTs, which we hope applies also to the CFT of the boson Hubbard model with reasonable 
accuracy. This theory depends upon two-dimensionless parameters: e and 7. However, both 
parameters can, in principle, be fixed by matching to the short-time, or u — > 00, limit of 
the correlators of the CFT of interest; e is related to Eoo, while 7 connects to a 3-point 
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correlator of the current J M and the stress-energy tensor [35J. 

Let us summarize our proposal for applying gauge-gravity duality to the CFT realized 
in the boson Hubbard model. We use effective field theory to motivate the effective action 



of the D + l dimensional gravity theory in Eq. (2.16). Conventional field-theoretic methods 
can be used to compute characteristics of the CFT at T = 0, and these can be matched to 
the gravitational theory to fix the values of parameters. Then we extrapolate to the long 
time limit at T > using the solution of the gravity equations of motion in the black brane 



background, as describe below Eq. (2.12). It is interesting that this procedure parallels 
that used for Boltzmann-like equations: use properties of the quasiparticles at T = to 
compute their scattering cross-section, so determine the collision term, and then solve the 
Boltzmann equation to extrapolate to the long time limit at T > 0. Evidently, the -D + l 
dimensional gravity theory has taken the place of the Boltzmann equation for CFTs without 
quasiparticle excitations. 

We are now ready to describe the results of this method in the computation of a(ou). The 
results depend upon the values of e and 7. However, the dependence on e can be scaled out 
by examining the ratio c(o;)/o' 00 . Furthermore, stability of the holographic theory requires 
that I7I < 1/12 [35]. Explicit results for this range of 7 are shown in Fig. [7j We see from 
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FIG. 7: Predictions [35] for the frequency-dependent conductivity cr(u) of a CFT in 2+1 dimensions 



from the holographic gravity theory in Eq. (2.16). The results shown depend only upon the 
parameter 7 which is restricted by stability requirements to the range I7I < 1/12. The 7 = 
case has a particle- vortex self-duality and applies to special supersymmetric models like the ABJM 
theory. 



Fig. [7] that the change in a(u) from the 7 = result to the maximal allowed values of 7 is 
quite limited: this stability is very encouraging, and is a partial a posteriori justification for 
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the validity of the gradient expansion applied to the holographic theory. 

A remarkable feature of Fig. [7] is that the results correspond very neatly to the conduc- 
tivities surmised from the Boltzmann equation in Figs. [4] and [5] the 7 > case is similar to 
the particle Boltzmann sketch in Fig. [4], while the 7 < case is similar to the vortex case in 
Fig. |5J Thus it is the sign of 7 which determines whether a given CFT is more accurately 
described by the particle-like or vortex-like excitations. We do not yet know the value of 7 
for the CFT of the boson Hubbard model, but the route is now open to its determination 
by the strategy we have outlined above. 

We close by noting a curious feature of Fig. [7J The 7 = case yields a frequency- 
independent a(u>). Recall that 7 = corresponded to the supersymmetric ABJM-like the- 
ories. These special theories have a sei/-duality under the particle-vortex transformation 
which accounts for this novel feature [36J. 

III. COMPRESSIBLE QUANTUM MATTER 
A. Condensed matter analysis 

We begin by defining compressible states, in a form applicable to both condensed matter 
and holographic theories [37]. 

• Consider a continuum, translationally-invariant quantum system with a globally con- 
served charge Q i.e. Q commutes with the Hamiltonian H. Couple the Hamiltonian 
to a chemical potential, /x, which is conjugate to Q: so the Hamiltonian changes to 
H — [iQ. The ground state of this modified Hamiltonian is compressible if (Q) changes 
smoothly function of /i, with d(Q)/dfi non-zero. 

Note that our definition of compressibility refers to the ground state, and so we need T = 0. 
Also, as noted in Section [TJ we will restrict our attention to compressible states in dimensions 
greater than unity (D > 2). 

Remarkably, it turns out there are only a few known states which satisfy these seemingly 
innocuous requirements. The states are: 

1. Solids: Translational symmetry is broken, and the matter 'crystallizes' into a periodic 
arrangement. Changing fi changes the period of the lattice, allowing a continuous 
variation in (Q). 

2. Superfluids: Here the global U(l) symmetry associated with conservation of Q is spon- 
taneously broken. These are gapless states which readily accept new particles into the 
Bose condensate, allowing a smooth variation in the density (Q). 

3. Fermi liquids: The simplest example of a Fermi liquid is the non- interacting Fermi 
gas. The fermionic particles occupy momentum eigenstates, with momenta inside a 
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Fermi surface which separates the occupied and empty states. This picture generalizes 
to interacting particles, to all orders in perturbation theory, as shown elegantly by 
Landau. In a Landau Fermi liquid, there is a sharp Fermi surface, and the only low 
energy excitations are sharp quasiparticles (renormalized from the bare fermions) in 
the vicinity of the Fermi surface. These quasiparticles carry opposite Q charges on 
either side of the Fermi surface. Changing fi changes the location of the Fermi surfaces, 
and the newly occupied states allow a smooth variation in (Q). 

Apart from mild variations which combine features of such states, these are the only examples 
of compressible quantum states in traditional condensed matter physics. 

The past decades have seen intensive search for other 'exotic' compressible states of 
condensed matter. As reviewed more completely in Ref. [201 this search has been motivated 
by the ubiquitous appearance of 'strange metal' behavior in a variety of correlated electron 
compounds, including the high temperature superconductors. From such studies a new 
class of compressible states appears to have emerged, which we refer to here generically as 
'non- Fermi liquids'. Many recent studies [2JJ 138H13] have emphasized the strong-coupling 
nature of the theory of non- Fermi liquids in two spatial dimensions (D = 3), and noted that 
important questions remain unresolved. 

Let us describe one of the simplest examples of a compressible non-Fermi liquid state. 
This was developed as a theory of a possible spin liquid (SL) state in frustrated quantum 
antiferromagnets [SHIE], and a complete derivation appears in the Supplementary Material. 
The primary degrees of freedom are fermionic 'spinons' f a , with a =t,4- a spin label. These 
are strongly coupled to an 'emergent' U(l) gauge field B^. Note that this gauge field is not to 
be confused with Maxwell gauge field associated with electromagnetism, whose fluctuations 
are ignored in all our analyses. It should also not be confused with the gauge field A a in 



Section II B which resides in the extended (D + l)-dimensional space. Rather, is a 
degree of freedom which emerges from the dynamics of the antiferromagnet, and encodes 
the complex quantum entanglement between valence bonds in the spin liquid. The action 
for this non-Fermi liquid state is simple (see Supplementary Material) 

d „ 1 



<S>sl = I d rdt 



fl[i^+B t -e F - ^-(V r - iB r ) 2 ) f a - B t J\f 



(3.1) 



where both f a and B^ are fluctuating functions of r and t, and the remaining parameters 
are constants. The Fermi energy ep controls the value of the fermion density (f a fa), but 
we are not free to choose its value. This density is coupled to a fluctuating U(l) gauge 
field which mediates long-range 'Coulomb' interactions, and so just as in the jellium model 
of the electron gas, stability of the system requires net neutrality in the gauge charge. We 
have thus included a background neutralizing charge density Af, and the value of Sf must 
be chosen so that 

{flf a )=M. (3.2) 
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Thus the system is not compressible with respect to the charge f^f Q ; moreover this charge 
is not associated with a global conservation law, but with a gauge invariance. Instead, the 
quantity realizing the globally conserved charge Q associated with compressibility is the spin 
density. The spin density has 3 components, and let us arbitrarily choose the z component, 
and so 

Q = fKpfa (3-3) 

where a z is a Pauli matrix. The "chemical potential" coupling to Q is the Zeeman coupling 
of an applied magnetic field. The compressibility of the SL state implies that its magneti- 
zation can be varied smoothly as a function of the applied magnetic field, and the magnetic 
susceptibility is non-zero at T = [45J. 

Like Landau Fermi liquids, a non-Fermi liquid SL state also has Fermi surfaces, but the 
character of the fermionic excitations near the Fermi surface is quite different. Formally, 
for interacting electron systems, the Fermi surface is defined by the momentum k = kp for 
which there is a zero in the inverse fermion Greens function: 



Gj\k = k F ,u = 0) = (3.4) 

Landau's Fermi liquid theory also predicts a specific singularity in the fermion Green's 
function near the Fermi surface: this singularity reflects the fact that the inverse-lifetime 
of the quasiparticles vanishes as the square of their distance from the Fermi surface. The 
energy of the quasiparticles is parametrically larger: it is linear in the distance from the 
Fermi surface, and so the quasiparticles are well-defined excitations. 

In the known non-Fermi liquids, the quasiparticles do not remain well-defined excitations 
away from the Fermi surface: they are strongly scattered by the low energy modes the U(l) 
gauge field. This is reflected in the singularity of the fermion Green's function near the 
Fermi surface, which has been argued to have the generic form [40] 



Gj\k,u) = q 1 -' 1 F(co/q z/2 ) (3.5) 

where we are expanding in the vicinity of the Fermi surface point at (kp, 0), with q x = k x —kp, 
q y = k y , and q — q x + q% (see Fig. j8|). The universal scaling function F is a complex- 
valued scaling function which determines the spectral density of the fermionic excitations. 
The parameters r] and z are critical exponents whose values can be estimated in various 
expansion methods [401 I4T] . However, the theory of these fermionic and gauge modes is 
strongly coupled [391 HQ] (in D = 3), and so the reliability of these estimates is unknown. 

We now make a few more remarks about the Fermi surfaces encountered so far: 
(i) The Fermi surface of non-Fermi liquids is sharply defined, even though the fermionic 



quasiparticles are not: the value of kp is precisely defined by Eq. (3.4), and the singularity 



of the scaling function in Eq. (3.5) on the Fermi surface 



[ii) The above results are expected to apply also to fermions coupled to non-Abelian gauge 
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+ q x 



FIG. 8: Fermi surface of a non-Fermi liquid. The "fuzzy" shading is used to indicate the broad 
spectral functions of the fermionic excitations. Note that despite this fuzziness, the position of the 
Fermi surface in momentum space is sharply defined. 

fields, or to order parameter fluctuations near a symmetry-breaking quantum phase transi- 
tion. Unlike the familiar case at zero density, there does not appear to be any fundamental 
difference between Abelian and non-Abelian gauge fields. 

(Hi) There is a crucial relation, known as the Luttinger relation, which constrains the total 
area enclosed by Fermi surfaces to the values of conserved charges. Specifically, let us con- 
sider a compressible state with Fermi surfaces labeled by I: each Fermi surface is associated 



with the singularity (3.4) in the Green's functions of fermions carrying global charge qt under 



symmetry associated with Q, and encloses area At. Then the Luttinger relation is (1H1 [50] 

J2qtA e = 4n 2 (Q). (3.6) 

i 

The fermionic excitations near the Fermi surface can also carry charges associated with 



fluctuating gauge fields, as in Eq. (3.1 ). If none of the Fermi surfaces have gauge charges, we 
usually obtain a Fermi liquid. If we have Fermi surfaces with and without gauge charges co- 
existing, we obtain a compressible state labeled as a fractionalized Fermi liquid [5U |52] . 



A Luttinger relation like Eq. (3.6) applies for each conserved U(l) charge Q associated 



with a compressible state. We also note [37] that for the case of an Abelian gauge field (as 



in Eq. (3.1)), there is an additional requirement for global neutrality in the gauge charge: 



in this case, there is additional Luttinger relation for the Fermi surfaces carrying the gauge 



charges, equating the sum of their areas to the background charge (which is M in Eqs. (3.1 ) 



and (3.2)). 



Our discussion of compressible states of condensed matter now suggests a conjecture [37] : 
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• All compressible states which preserve translational and global U(l) symmetries must 
have Fermi surfaces, but they are not necessarily Fermi liquids. The areas enclosed by 
these Fermi surfaces obey a Luttinger relation for each conserved charge. 

Note that we have not imposed any requirement that our underlying Hamiltonian have any 
fermionic degrees of freedom. It could be defined solely in terms of bosons, or have both 
fermions and bosons; recent examples of boson-only models are in Refs. EE1 E3J EH There 
are fermionic excitations near each Fermi surface, and these fermions could be fractions or 
composites of the underlying particles. 

The rationale behind this conjecture is that (in D — 1 spatial dimensions) we need a 
surface of dimension D — 2 of zero energy excitations to have adequate phase space to allow 
(Q) to vary smoothly as a function of the chemical potential. In D > 2, bosons can't 
generically have such a surface of zero energy excitations: there will be negative energy 
states on one side of this surface, rendering the system unstable to Bose condensation. In 



contrast, fermions are allowed to have such surfaces with a singularity as in Eq. (3.4), as is 
evident already from the free fermion theory: the negative energies now represent hole-like 
excitations of the fermions. 

We are finally in a position to state the challenge to the holographic approach to con- 
densed matter physics. Can gauge-gravity duality provide a classification of possible states 
of compressible quantum matter? Clearly, any such classifications should include the famil- 
iar phases: solids, superfluids, and Fermi liquids. Are there additional states, and do they 
correspond to the non-Fermi liquid states with Fermi surfaces described above? If so, we 
can hope that they will provide a new perspective on their strong-coupling physics. Or is 
the conjecture incorrect, and are there entirely new types of compressible quantum states 
which have been overlooked in condensed matter studies? 



B. Holographic analysis 

The study of compressible states using holography is the focus of much current research 
[T0HT8] . This research does not yet have definitive answers to the questions posed at the 
end of the previous section. However, rapid progress has been made recently, and here I will 
give my perspective on the current state of the theory. This research promises to eventually 
achieve a holographic classification of the possible states of compressible quantum matter. 
Such a classification will be useful for condensed matter applications, especially in two spatial 
dimensions, regardless of the artificial nature of the microscopic degrees of freedom used to 
realize the compressible phases. 



The basic starting point is to extend the Einstein-Maxwell theory Sem in Eq. (2.8) which 



we used in the conformal case to a situation with non-zero density. Given the relationship 



in Eq. (2.12), we can turn on a finite density by including a chemical potential, in which 
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case we now have the boundary condition 



AJr, t, u — y oo) 



(3.7) 



while the other components of the gauge field vanish in the corresponding limit near the 
boundary. We will now work to all orders in /i, rather than the linear response in the source 
term in Section III Bl 

In the simplest (and frequently used) approach, we make no further change. We solve 



the equations of motion of the action in Eq. (2.8) subject to the boundary condition in 



Eq. (3.7) for the gauge field. Such a solution is of the Reissner-Nordstrom form [55]: the 



metric remains as in Eq. (2.13), but the function f(u) in Eq. (2.14) is replaced by 



f(u) = 1 + 



L 2 kV 
2e 2 R 2 



- 1 + 



L 2 K 2 fi 2 \ fR 
2e 2 R 2 ) [u 



(3.8) 



Again R is the position of the horizon, as it can be checked that f(R) = 0. In addition, we 
have a non-zero gauge field given by 



A t (u) =n[l- 



R 



u 



(3.9) 



which obeys Eq. (3.7), and vanishes on the horizon. As above Eq. (2.15), we can compute 



the temperature of the horizon, and now find that it equals 



, m 3hvR ( L 2 h 2 k 2 
k B T = , ^ I 1 



4ttL 2 



6e 2 R 2 



(3.10) 



as in Eq. (2.15), this equation fixes R in terms of physical parameters. A sketch of this 



geometry, generalizing Fig. |6j is shown in Fig. [9j Note that the total 'electric' flux in the 
bulk is conserved between the horizon and the boundary at u = oo: this follows from the 
potential in Eq. (3.9), which shows that the electric field ~ 1/u 2 , and the metric in Eq. (2.13), 
which shows that the surface area ~ u 2 . 

This gravitational solution can now be applied to examine physical properties of the 
quantum system on the boundary as a function of the temperature, T, and the chemical 
potential, \i. For T ^> \i this was examined in Ref. [SBJ and successfully used to develop a 
general theory of thermoelectric transport in compressible quantum systems near quantum 
critical points. We will not review this here, as our primary interest is in the limit T <C \x for 
non-superfluid states. I also note in passing the interesting recent work of Bhattacharya et al. 
[5T] which has obtained new insights on the low temperature hydrodynamics of superfluids 
using the dual gravity approach. 

A notable feature of the T — > limit of the above Reissner-Nordstrom solution is apparent 
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Near-horizon 
metric is 
AdS, x R 2 



FIG. 9: Schematic (adapted from Ref. ITT]) of the proposed geometry of gauge-gravity duality for 
compressible quantum matter, generalizing Fig. [6] to a AdS4-Reissner-Nordstr6m black brane. The 
non-zero density (Q) on the boundary sources an electric flux which extends into the bulk to the 
horizon. Consequently, there must be compensating charges sourcing the electric flux 'behind' the 
horizon. 



from Eq. (3.10): the horizon radius, R remains finite when /x is non-zero, with 



R 



\7Et 



-/' 



T = 0. 



(3.11) 



This is different from the conformal case of Section II B , where the horizon was absent at 
T 



0. An immediate consequence is that our quantum system on the boundary, has a 
non-zero entropy density in its ground state, as can be verified by a computation of the free 
energy of the gravitational solution. The free energy computation [12] also shows that the 
compressibility of the ground state is finite, where the charge density is a smooth function 
of /i: 



(Q) 



T = 0. 



(3.12) 



VEePL 

So we have succeeded in obtaining a phase of compressible quantum matter, but it has 
the unacceptable feature of a ground state degeneracy which increases exponentially with 
system size. 

However, we will not discard this solution right away: let us examine some of its properties 
more carefully and see if a physical interpretation emerges. This will also help point the 
way towards improving this holographic description. 



We examine the form of the metric more closely; using the relationship (3.11 ) in Eq. (3.8) 
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we find 

,/ ^ ( R\ 2 ( 2 R 3i? 2 \ 

/(„= !__ ! + _+ ; T = 3.13 

\ u J \ u u z J 

Note that f(u) has a double zero at the horizon u — R, unlike the conventional single zero 
in the conformal case in Eq. (2.14). Such a horizon is known as extremal, and this feature 
is linked to its anomalous properties. Defining a co-ordinate which is zero at the horizon, 
u = u — R, and expanding to the lowest non- vanishing order in u, the metric in Eq. (2.13) 
becomes 

* , -Ks)'*-® , « ,+ (i) v 

The notable feature of this metric is its separation into the first two terms dependent only 
upon t and u, and the last term dependent only upon the spatial co-ordinate. This means 
that spacetime has factorized into a spatial R 2 , and a curved space in u and t. By comparison 
with Eq. (2.9), the reader will recognize that the latter space is AdS 2 , and so we conclude 
that the near-horizon spacetime of a Reissner- Nordstrom black brane is AdS 2 x R 2 . This 
factorized form of the metric has strong consequences for all the low energy properties of 
the compressible quantum system on the boundary. 

Let us begin by interpreting the AdS 2 factor. What kind of quantum system does this 
describe? As I have reviewed in more detail elsewhere [TOJ [TTj, [58] , it describes the physics of 
a single quantum impurity universally coupled to a CFT, analogous to the overscreened fixed 
points of multichannel Kondo problems jSH], or a vacancy in a two-dimensional insulating 
antiferromagnet at a magnetic ordering critical point [60]. There is a conformal structure 
to the correlations near the impurity, and the ground state has a finite entropy in the limit 
T — > [SHJ EI]. All these features are physical and robust, and potentially realizable in 
generic experimental systems. Supersymmetric generalizations of such quantum impurity 
systems can be solved [62] by gauge-gravity duality, and explicitly reveal an AdS 2 metric in 
the low energy physics. 

We are now ready to interpret AdS 2 x R 2 . The factorized form of the metric means that 
we can view the two-dimensional quantum system as consisting of an infinite number of 
quantum impurity problems, one at each position in space. Because there is an ultraviolet 
cutoff above which the AdS 2 factorization does not hold, there is a similar cutoff in spatial 
separation below which the quantum impurities remain coupled. Thus the low energy physics 
is controlled by a finite density of independent quantum impurities. Each quantum impurity 
is coupled to a gapless environment, which is presumably a mean-field description of the 
other quantum impurities. Many condensed matter theorists will immediately recognize the 
similarity of this picture to that of dynamical mean field theory (DMFT) [63] • Such large 
dimension/co-ordination number approximations have also been applied to the vicinity of 
magnetic ordering transitions, and descriptions of compressible non- Fermi liquid states have 
been obtained [EU I64TI68"] . Such states now have a non-zero ground state entropy density 
a consequence of adding up the entropy of the finite density of decoupled quantum 
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impurity problems. It was argued [58] that these compressible critical states, obtained in 
the limit of large spatial dimension or lattice co-ordination number, which provide a specific 
microscopic realization of the physics of AdS2 x R 2 . The holographic theory shares the 
feature of having a non-zero ground state entropy density, and we will now see that the 
structure of correlations of local operators also match. 

Correlations of the boundary theory are determined by inserting probe fields in the 'bulk' 
gravity theory, and computing the limiting values of their Green's functions near the bound- 
ary. As representative examples of such correlations, we introduce the simplest example of 
'matter' fields in the gravity theory: a complex scalar field 0, and a Dirac fermion, The 
action of these fields respectively has the schematic form 



S b = J d A x^~g [-|V</» - iAcf)\ 2 - m 2 \<f)\ 2 ] 



S 



1 

— c 
4 



r ■ ( d + ^coabT^ -iA)m - mmm 



(3.15) 



where oo a b is the spin connection of the curved geometry, and the T represent Dirac matrices. 
The parameter m determines the scaling dimension of the operator being used to probe 
the compressible quantum states. We are now faced with the problem of determining the 



inverses of the differential operators in Eq. (3.15) in the 4-dimensional space defined by 



Eqs. (2.13) and (3.13), and taking their limit u — > oo limit. This is an intricate problem 

In 



in solving differential equations on a curved space, addressed in recent work 
the low frequency (oo) limit, the Green's functions of the boundary compressible quantum 
theory had the following form for both the boson and fermion cases [72J: 



G-\k,oo) =C(k) + D(k)u Uk 



0. 



(3.16) 



where C(k), D(k), and are smooth functions of the wavevector k. This peculiar behavior, 
with a smooth dependence on k, and a singular dependence on oo, is a direct consequence 



of the AdS2 x R factorization of the geometry in Eq. (3.14): the independent quantum im- 



purities have temporal correlations which decay with a power-law in time, but are spatially 



uncorrelated. The behavior in Eq. (3.16) does not correspond to any of the generic com- 



pressible condensed matter phases considered in Section |III A| Instead, it is essentially the 
behavior observed in the compressible critical states in the limit of large spatial dimension or 
lattice co-ordination number [61, 64-68] . The correspondence between these large dimension 



models and the holographic solutions extends also to the extension of Eq. (3.16) to non-zero 
temperatures. 



The computation of the inverses of the operators in Eq. (3.15), also yields stability con- 



ditions on the gravity theory. For the boson action, Sb, the condition is C(k) > 0. This 
is satisfied for m sufficiently large. However, with decreasing m there is instability towards 
condensation of the scalar, leading to a new compressible phase with superfluidity 
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we will discuss this phase further in Section HI B 1 



For the fermion case, C(k) is allowed to be either positive or negative. If C(k) changes 
sign as a function of the k, we obtain a Fermi surface [72] at k — kp, where C(kp) = 0, 
by Eq. (3.4). The singularity of the Green's function in Eq. (3.16) at the Fermi surface is 
a special case of the generic form in Eq. (3.5); the latter is singular as a function of both cu 
and k — kp, and reduces to the singular part of the present form obtained for the AdS2 x R 2 
theory only in the limit where the dynamic critical exponent z — > oo. Indeed, the value 
z = oo is another way of characterizing the peculiar phase of the holographic theory with a 
geometry with factorization between space and time. 

Finally, we should compare the compressible state described by the Reissner-Nordstrom 
geometry to the conjecture towards the end of Section HI A We have obtained a compressible 
state which preserves both the translational and global U(l) symmetry associated with Q. 
We found Fermi surfaces only over a limited parameter regime, but even when present, they 
enclose areas which do not obey the Luttinger relation in Eq. (3.6). For the conjecture to 
be valid, the deficit in the Luttinger relation must be made up by 'hidden' Fermi surfaces 
associated with fermions which carry additional gauge charges, as in the non-Fermi liquid 
and fractionalized Fermi liquid states discussed in Section HI A The correlators of fermions 
with gauge charges are not readily computable in the gauge gravity duality, and are a worthy 
topic for future research. 



1. Beyond AdS 2 x R 2 

Given the non-zero ground state entropy density of the compressible state with z = oo 



theory described in Section III B , it seems clear we need to move beyond the Reissner- 



Nordstrom geometry obtained via the equations of motion of the Einstein-Maxwell action in 



Eq. (2.8). In particular, at the very least, we should include the back-reaction of the matter 



fields in Eq. (3.15) upon the geometry. This is an active topic of current research, and I 
summarize a few recent developments. 

As an example of a back-reaction on the metric, we can move to the superfluid state 
where the bosons are condensed, and recompute [801 EE] the metric g and gauge field A t 
from the classical equations of motion of Sem + <Sb, where is treated as a c-number as in 
the Bose-Einstein theory at T = 0. In the non-superfluid state, some studies [8"2~H8"6] have 
accounted for the back reaction of fermions in the Thomas-Fermi approximation, where the 
fermion density at a given u is a function only of the local chemical potential. In both of 



these qualitatively similar back-reaction is obtained, which is illustrated in Fig. [L0 

The electric field is screened by the bosonic or fermionic matter in the bulk, so that the 
electric field vanishes as u — > 0. In the opposite limit u — > oo, the electric field remains 
pinned by the value of (Q), as is required by Gauss' Law on the boundary. For the metric, 
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IR Lifshitz 
metric 



FIG. 10: Schematic of the geometry for compressible quantum matter, after accounting for back- 



reaction from the bulk matter in Eq. (3.15) in a Thomas- Fermi or Bose-Einstein theory. There 



is no horizon, and the small u geometry has a 'Lifshitz' form, as in Eq. (3.17). The bulk charge 
density screens the electric field, so that the electric field vanishes as u — > 0. 



the following qualitative structure is obtained as u — > 



ds 2 



du 2 



a- 



- u 2z dt z + u 2 dr< 



(3.17) 



where a and z are constants. It has become the practice in the string theory literature to 
call this the 'Lifshitz' metric [SZ] (for not particularly good reasons). A key property of this 



metric is that it is invariant under the following transformations, which generalize (2.3) and 



(2.10) 



t^t/b z , r^r/b , u^bu, 



(3.18) 



and identify z as the dynamic critical exponent. The AdS4 metric of Eq. (2.9) corresponds 
to z = 1. The AdS2 x R 2 metric of Eq. (3.14) is obtained in the limit z — > oo (as expected), 
after setting u = u 1 ^. 

The existence of this scale-invariant structure in the small u limit is puzzling, and is 
not well understood. It implies there are emergent low energy excitations which scale with 
the dynamic critical exponent z. In the traditional condensed matter superfluid state, the 
only low energy excitations of a superfluid are the Goldstone modes, and there are no 
dynamic critical fluctuations. This indicates that the holographic superfruids [HJ[75] are n °t 
conventional superfluids at T = 0, and perhaps co-exist with a non- Fermi liquid sector. For 
the non-superfluid compressible state obtained using the Thomas-Fermi theory, it is possible 
that these critical excitations are associated with the gapless gauge excitations which led to 



Eq. (3.5) in the non- Fermi liquid. However, the latter excitations are associated with long 
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wavelength, low energy fluctuations at momenta near the Fermi surface, and these are not 
contained in the Thomas-Fermi theory. Indications from recent work [88H92] are that these 
low-energy excitations arise from the presence of an infinite number of Fermi surfaces, with 
a Fermi wavevector which varies continuously with u. 

For condensed matter applications, we therefore need to move the theory of compressible 
states to a regime with a small number of Fermi surfaces. This was addressed recently 
in Ref. [92] by filling up fermionic states of <S/ in a suitable metric. A key point made in 
Ref. 92 is that a Fermi liquid must be a confining state of the gauge theory defining the 
CFT, as we discuss in the Supplementary Material. In holographic studies of zero density 
theories of particle physics, it is known that a confining state corresponds to a geometry 
which terminates at finite u = u m , rather than extending all the way to u = 0; see Fig. [TT 



Assuming such a geometry, Ref. [92] discussed a complete quantum solution of the theory 



IR confining 
metric 




FIG. 11: Confining geometry for a Fermi liquid state. The geometry terminates at u = u m and 
the bulk charge density screens the electric field, so that the electric field vanishes as u — > u m . 



Sf in Eq. (3.15), without making the Thomas-Fermi approximation. A conventional Fermi 
liquid state was found, without extraneous excitations. Moreover, the Luttinger relation in 



Eq. (3.6) was found to be exactly obeyed with a small number of Fermi surfaces, including 



the simplest case with a single Fermi surface; this Luttinger relation was a consequence 
of Gauss's Law for the electric flux in the u direction. Confining geometries were also 
considered for the superfluid state [93] US] , and our present arguments indicate that these 
are traditional superfluids. 

Holographic realizations of the familiar compressible phases of condensed matter physics, 
Fermi liquids and superfluids, thus finally appear to have been found. Extensions of this 
understanding to non- Fermi liquids and fractionalized Fermi liquids are important topics for 
future research. The connection between the Luttinger relation and Gauss' law [92] indicates 
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that the holographic realizations of these phases will have at least part of the electric flux 
'leaking out' at u = 0, as in Fig. |9j and in a recent study [89] . 

IV. CONCLUSIONS 

We divided our discussion into two classes of systems realizable in condensed matter 
physics: conformal and compressible. 

Conformal systems were discussed in Section |TTJ These concerned generic models which 
realize conformal field theories (CFTs) at a quantum phase transition, the simplest being the 
superfluid-insulator transition of bosons at integer filling in a periodic potential in two spatial 
dimensions. Many T = properties of the CFT are accurately computable by traditional 
field theoretic perturbative expansions. However, this success does not extend to long-time 
correlations at T > 0, because the orders of limit of t — > oo and T do not commute 
for all such methods. The Boltzmann equation, and its many descendants, are traditionally 
used to resum the perturbative methods, and access the long time limit; these are only 
expected to be reasonable when particle-like excitations are long-lived. The holographic 
method provides an alternative description of relaxational and transport properties at long 
time: it does not assume any picture of particle-like excitations and so is far removed from, 
and complementary to, the Boltzmann description. We described a gradient expansion 
method [35J for generating an effective holographic theory for generic CFTs. The coupling 
constants of this effective theory can be fixed by matching to various T = correlators of the 
CFT. Then the effective theory generates useful results for the long-time limit, including the 
values of transport co-efficients. I believe this method offers hope for quantitative, testable 



predictions not only for the linear transport problems discussed in Section |II B[ but also 
non-linear and non-equilibrium problems in the vicinity of quantum critical points. 
Compressible systems were studied in Section HI A| The familiar condensed matter ex- 



amples are superfluids, solids, and Fermi liquids (and their variants). An exotic class of 
compressible states have emerged in condensed matter studies motivated by the strange 
metal problem of correlated electron systems [20J. These are varieties of non- Fermi liquids, 
all of which have sharp Fermi surfaces, but the fermionic quasiparticles near some or all of 
the Fermi surfaces are not well-defined because of their strong coupling to deconfined gauge 
fields, or to gapless modes at a quantum phase transition. This led us to the conjecture [57] . 



at the end of Section HI A that all compressible, continuum quantum systems which do not 
break translational or the global U(l) symmetry (associated with the conserved density) at 
T = 0, must have Fermi surfaces; further, the areas enclosed by these Fermi surfaces must 



obey the Luttinger relation in Eq. (3.6) 



The holographic method has the potential to provide us with an alternative classification, 
and a deeper understanding of the compressible states of quantum matter. Such a classi- 
fication should surely include the familiar solids, superfluids, and Fermi liquids, and it is 



27 



not yet settled if the exotic compressible states obtained via holography correspond to those 
obtained so far in the condensed matter literature. It remains possible that the Fermi sur- 
face conjecture is false, and holography leads us to qualitatively new types of compressible 
phases. 

The Reissner-Nordstrom solution of the Einstein-Maxwell action provides the simplest 
theory of a compressible state. We reviewed the physical properties of this state: a non-zero 
ground state entropy, 'locally critical' correlations in time, a non-Fermi liquid Fermi surface 
with dynamic critical exponent z = oo, and T > correlations similar to those of critical 
quantum impurity problems. We argued that all these properties were in close correspon- 
dence with the 'fractionalized Fermi liquid' phase [5TL 152] of lattice quantum models solved 
in the limit of large spatial dimensionality or large lattice co-ordination number [58J . In Sec- 



tion HI B 1 we discussed ongoing research moving beyond the Reissner-Nordstrom geometry. 



In many studies, back-reaction of the matter fields on the gravity theory was examined in 



Thomas- Fermi-like approximations, and led to a low energy 'Lifshitz' geometry in Eq. (10) 
with a finite z. However, the physical interpretation of the critical low energy excitations 
with such a geometry remains unclear: such excitations are absent from conventional super- 
fluids and Fermi liquids, and it is unlikely they correspond to the modes of known non- Fermi 
liquids. They may well be artifacts of the presence of an infinite number of Fermi surfaces, 
with a continually varying Fermi wavevector, along the emergent holographic direction [92J. 
In other studies [92H91] . confining geometries which terminate along the holographic u di- 
rection have been considered, and these do realize traditional superfluids and Fermi liquids, 
without extraneous excitations. These studies also indicate that the compressible non- Fermi 
liquid (or fractionalized Fermi liquid) phases will be realized in geometries with all (or part) 
of the electric flux remaining unscreened until the end of the holographic direction. 

It does seem that further refinements of the theory are needed before the links between the 
condensed matter and holographic theories of compressible quantum matter are completely 
established. Given the pace of progress in the past few years, we can hope that many of 
these issues will be resolved in the not too distant future. 

Note added: A new approach to holographic theories of compressible quantum matter has 
since appeared in Ref. 191)1 
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Abstract 

Field theories of conformal and compressible states of matter 
are derived from lattice Hubbard models. 



These notes are adapted from earlier reviews [U E] • 



I. CONFORMAL QUANTUM MATTER 



We will demonstrate that the superfluid-insulator quantum phase transition of the boson 
Hubbard model in two spatial dimensions is described by a conformal field theory (CFT). 



The boson Hubbard model is (as in Eq. (2.1) in the main paper) 



(ij) i i 

where hi is the canonical boson annihilation operator, n% = b\bi is the boson number operator, 
w is the hopping matrix element between nearest-neighbor sites, U is the on-site repulsive 
energy between a pair of bosons, and \i is the chemical potential. Let us assume that the 
average boson density is exactly n per site, where n is a positive integer. For U/w 3> 1, 
the ground state is simply 

\GS) = H(b\r\0), (1.2) 

i 

where |0) is the empy state with no bosons. In the same limit, the lowest excited states are 
"particles" and "holes" with one extra or missing boson, 

\Pi) = b]\GS} } (1.3) 
\ht) = bi\GS). (1.4) 



1 



For w/U = strictly, the particle and hole energies (relative to the ground state) are 



£j°) = Un - fi, 4 0) = U(l - n ) + p, (1.5) 

For < w/U < 1, these states will develop dispersion. By considering the first order 
splitting of the degenerate manifold of particle or hole states (degeneracy associated with 
the site of the particle or hole), one obtains (considering the square lattice for simplicity) 

k 2 

E p (k) = - 2w(n + l)(cosfc x + cosfcj,) « A p + 



2m p 
k 2 

E h (k) = E ( h 0) - 2wn (cos k x + cos k y ) « A h + - — , (1.6) 

where we have Taylor expanded around the minimum at k — 0, giving m p = l/(2(no + l)w) 
and rrih = l/(2now) The excitation gaps, A p /h, are 

A p = Un — /j, — Aw(n + 1) 

A h = U(l-n )+p- 4wn . (1.7) 

to this order in w/U. As long as both of the these gaps are positive, our starting point of a 
Mott insulating state with an average of n particles per site is stable. 

When one of the gaps vanish, the Mott insulator is no longer stable, and we have a 
quantum transition to a superfluid state. Let us assume that it is A p that vanishes first 
with increasing w. The transition then corresponds to a Bose-Einstein condensation of 
particles, with — A p acting as the effective chemical potential. At T = 0, an increasing 
chemical potential implies an increasing particle density, and so the superfluid state will 
have a density greater than that of the Mott insulator. Similarly, if the value of \x is such 
that A h vanishes first, the superfluid state will have a density smaller than that of the Mott 
insulator. 

However, let us consider the special case when the density of both the Mott insulator and 
the superfluid are equal to no; this is often naturally the case under experimental conditions. 
Our reasoning makes it clear that this is only possible if \x is chosen so that A p = Ah = A. 
This both gaps vanish simultaneously, the insulator-superfluid transition corresponds to 
condensation of both particles and holes (which can be viewed as "anti-particles"). This 
symmetry between particles and anti-particles is responsible for the relativistic structure of 
the low energy theory. 

Let us now proceed to derive the effective action for the low energy theory near the 
insulator-superfluid transition. While it is possible to derive a field theory of this condensa- 
tion from Hb, we instead just write it down based on our simple physical picture. We model 
the particle and hole excitations by fields p(r, r), h(r, r) respectively, in the imaginary time 
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(r) path integral. The weight in the path integral is, as usual, the Euclidean action, 



S h 



drd r 



P 



d_ 



+ A 



2m r 



P 



d V 2 

dr 2mft 



J h - A(p j h j + ph) + 



Here we have included a term A which creates and annihilates particles and holes together in 
pairs, which is expected since this conserves boson number. Microscopically this term arises 
from the action of the hopping w on the naive ground state, which creates particle-hole pairs 
on neighboring sites, so A ~ 0(w) (the spatial dependence is unimportant for the states 
near k = 0). We have neglected - for brevity of presentation - to write a number of higher 
order terms involving four or more boson fields, representing interactions between particles 
and/or holes, and other boson number-conserving two-body and higher-body collisional 



processes. Note that the dependence upon w/U in Eq. (1.8) arises primarily through implicit 
dependence of A. 

Without loss of generality, we assume A > 0, and change variables to the linear combi- 
nations 



1 



V2 



;i.9) 



Then the quadratic terms in the action are 



Sb = J drd 2 r 
1 



e t^_ e ^! + (A _ A)H 2 +(A+A)|e|2 



dr 



dr 



4m p 4:171^ 



(|W| 2 + |V£| 2 ) + 



4m p 



. (1.10) 



Notice that the quadratic form for ip becomes unstable, before that of £. So let us integrate 
out £, expanding the resulting action in powers and gradients oiip. In this manner we obtain 
the theory for the superfluid order parameter ip [3] 



drd 2 r 



1 



(A + A) 



dijj 



+ (A-A)M 2 + 



1 



1 



4m„ 4m/j 



;i.ii) 



This is the promised relativistic field theory, written in Eq. (2.2) in the main paper. The 



energy gap for both particle and hole excitations is a/A 2 — A 2 , and this vanishes at the 
quantum critical point. 

There are a number of additional higher-order terms, not displayed above, which are not 
relativistically invariant. However, all of these are formally irrelevant at the Wilson-Fisher 
fixed point which controls the critical theory. 

Finally, we note that the scaling properties and relativistic invariance of the critical point 
are sufficient to establish its invariance under conformal transformations. 
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A. Quantum critical transport 



To illustrate the general issues, we begin by computing the transport properties of the 
free field theory of a complex scalar with mass m, written in a Lorentz invariant notation: 

S^ = J d D r[\d^\ 2 + m 2 \^\ 2 } (1.12) 

This theory can be obtained from Eq. ( |1.11[ ) at u = 0, after appropriate rescalings of co- 
ordinates and fields. 

The conserved electrical current is 



■1, 



;i.l3) 



Let us compute its two-point correlator, K^ik) at a spacetime momentum fc M at T = 0. 
This is given by the one-loop diagram which evaluates to 



KUk) 



(Pp / (2 Pll + k v )(2p v + k v ) _ 



^ 3 V {{p + k) 2 + m 2 )(p 2 + m 2 ) p 2 + m 2 



1 ( ^ k^k v 



in 



'fJ,U 



dx- 



k 2 {l-2xf 



;i.i4) 



k 2 J Jo ^/m 2 + k 2 x(l - x)' 
The second term in the first equation arises from a 'tadpole' contribution which is omitted 



in Eq. (1.13). Note that the current correlation is purely transverse, and this follows from 



the requirement of current conservation 



k K 



0. 



;i.i5) 



Of particular interest to us is the K 00 component, after analytic continuation to Minkowski 
space where the spacetime momentum k^ is replaced by (u, k). The conductivity is obtained 
from this correlator via the Kubo formula 



a(uj) = lim — 
k^o k 2 



K 00 (uj,k). 



(1.16) 



In the insulator, where m > 0, analysis of the integrand in Eq. (1.14) shows that that 
the spectral weight of the density correlator has a gap of 2m at k = 0, and the conductivity 



in Eq. (1.16) vanishes. These properties are as expected in any insulator. 

At the critical point, where m = 0, the fermionic spectrum is gap less, and so is that of 



the charge correlator. The density correlator in Eq. (1.14) and the conductivity in Eq. (1.16) 
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evaluate to the simple universal results 



K 00 (u, k) 



16 
1 

16' 



;i.i7) 



Going beyond the free field theory in Eq. (1.12), the effect of interactions can be accounted 



for order-by-order in u. In the renormalization group approach, u takes the value specified 
by the Wilson-Fisher fixed point at the quantum critical point. Combined with the absence 



of of divergencies in the perturbative expansion (which is a consequence of Eq. (1.15)), this 



means the only effect of interactions is to change the pre-factor in Eq. (1.17) to a different 
universal numerical value. So we write 



K Q0 (u), k) = a a 
a(uS) = a 



Vk 2 - tu 2 



(1.18) 



where is a universal number dependent only upon the universality class of the quantum 
critical point, whose value can be computed by various expansion methods. 



1. Non-zero temperatures 



We begin by repeating the above computation for the free field theory at T > 0. This 



only requires replacing the integral over the loop frequency in Eq. (1.14), by a summation 



over the Matsubara frequencies which are quantized by integer multiples of 2irT. Such a 



computation, via Eq. (1.16) leads to the conductivity 



He[<7(w)] 
V 



VS(u) 



2m) 



16 

k 3 dk 



Am' 



4a; 2 
1 



coth 



\UJ\ 

4T 



8T J (k 2 + m 2 ) sinh 2 (Vk 2 + m 2 /2T) 



1.19) 



the imaginary part of <j(uj) is the Hilbert transform of Re[cr(co>)] — 1/16. Note that this reduces 



to Eq. (1.17) in the limit u> T. However, the most important new feature of Eq. (1.19) 



arises for u ^ T, where we find a delta function at zero frequency in the real part. Thus the 
d.c. conductivity is infinite at this order, arising from the collisionless transport of thermally 
excited carriers. This is clearly an artifact of the free field theory. 



At non-zero u, collisions between carriers invalidate the form in Eq. (1.19) for the density 



correlation function, and we instead expect the form dictated by the hydrodynamic diffusion 
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of charge. Thus for K 00 , Eq. (1.18) applies only for uj ^> T, while 

Dk 2 



K 00 (uj,k) = x 



Dk 2 — itu 



uj <T. 



;i.2o) 



Here \ is the charge susceptibility (here it is the compressibility), and D is the charge 
diffusion constant. By the universality of the Wilson-Fisher fixed point, we expect that 
these have universal values in the quantum critical region: 



X = C X T 



D 



C 



D 



T ' 



1.21) 



where again C x and Co are universal numbers. For the conductivity, we expect a crossover 
from the collisionless critical dynamics at frequencies w > T, to a hydrodynamic collision- 
dominated form for w <T. This entire crossover is universal, and is described by a universal 
crossover function 

a{u) = E{uj/T). (1.22) 



The result in Eq. (1.18) applies for oj 3> T, and so 



E(OO) = (Too. 



;i.23) 



For the hydrodynamic transport, we apply the Kubo formula in Eq. (1.16) to Eq. (1.20) and 
obtain 

E(0) = C X C D (1.24) 



which is a version of Einstein's relation for Brownian motion. 



II. COMPRESSIBLE QUANTUM MATTER 

Now we will consider the Hubbard model for fermionic particles with spin S = 1/2 
(electrons) on the triangular lattice. 

For small U/w, the ground state of this model is a metal, rather than a superfluid. This 
is because the fermions cannot condense; instead they occupy all single particle states inside 
a 'Fermi surface' in momentum space, forming a Fermi liquid. 

For large U/w, and with a density of one electron per site, we do expect an insulating state 
to form, with a gap to both particle and hole excitations, just as was the case for bosons. 
However, the electron localized on each site of the lattice now has a spin degeneracy, and we 
also have to specify the spin wavefunction in the insulator. At the largest values of U/w, it 
is believed that the insulator has long-range antiferromagnetic order; we will not study this 
ordered state here. The nature of the insulator at smaller U/w, and in particular, in the 
vicinity of the insulator-metal transition is still a question of some debate. In the following, 
we will assume that the insulating state proximate to the critical point is a particular "U(l) 
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spin liquid" , which we will describe more completely below. 
The Hubbard Hamiltonian is 



H — "' ^ \ /n' l" 



E 



-fi (ntf + Tin) + U [ntf 



(2.1) 



Here Cj a , a =j, 4- are annihilation operators on the site i of a triangular lattice. The density 
of electrons is controlled by the chemical potential \x which couples to the total electron 
density, with 

(2.2) 



For completeness, we also note the algebra of the fermion operators: 



0. 



(2.3) 



Let us begin by considering the case U — 0. Then the ground state is a metal at all 
densities, with a Fermi surface separating occupied and empty states in momentum space. 
Landau's Fermi liquid (FL) theory describes how the free-electron model of a metal can be 
extended to non-zero U. For our purposes, we need only two basic facts: (i) the fermionic 
excitations near the Fermi surface are essentially non-interacting electrons, and (ii) the area 
enclosed by the Fermi surface is equal to the electron density — this is Luttinger's theorem, 
which we state more explicitly below. 

At U = 0, the Hamiltonian of the FL metal is 



H = ^ c^(fc) — fx — 2t^cos(fc ■ e\) + cos(fc ■ e 2 ) + cos(fc • e 3 ) \ 



cJk) 



(2.4) 



where the are as shown in Fig. [T] The reciprocal lattice consists of the vectors G 
Yji n iGii where 



An 

Gi = —{e x - e 2 ) 



An 

G 2 = y(e 2 - e 3 ) 



47T 

^3 = y(e 3 - e Xl 



(2.5) 



The electronic dispersion in Eq. (2.4) is plotted in Fig. 2 it only has simple parabolic minima 



at k — 0, and its periodic images at k = G, and there are no Dirac points. At any chemical 
potential, the negative energy states are occupied, leading to a Fermi surface bounding the 
set of occupied states, as shown in Fig. [3j Luttinger theorem states that the total area of 
the occupied states, the shaded region of the first Brillouin zone in Fig. [3] occupies an area, 
A, given by 

A/", (2.6) 



A_ 

2vr 2 



where Af = J2 a <^ a c a is the total electron density. This relationship is obviously true for 
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FIG. 1: The triangular lattice 




FIG. 2: The electronic dispersion in Eq. (2.4) for fi = and t = 1 



free electrons simply by counting occupied states, but it also remains true for interacting 
electrons. 

Now we turn up the strength of the interactions, U, at a density of one electron per 
site. By the same argument as that for bosons, an insulator will appear for sufficiently large 
U . As stated above, we will focus particular route to the destruction of the small U Fermi 
liquid, one which reaches directly to an insulator which is a 'spin liquid' [5H7j. The spin 
liquid insulator is a phase in which the spin rotation symmetry is preserved, and there is a 
gap to all charged excitations. However, there are gapless spin excitations, and an emergent 
compact U(l) gauge field in a deconfined phase. 

The key to the description of this metal-insulator transition is an exact rewriting of the 



Hubbard model in Eq. (2.1) as a compact U(l) lattice gauge theory. To derive this, let us 



S 




FIG. 3: The Fermi surface of Eq. (2.4) for /x = 1/2 and t = 1; the occupied states are shaded. Also 



shown are the periodic images of the Fermi surface in their respective Brillouin zones. 

proceed with using the same strategy as that used in Section |T] for the boson Hubbard model. 
So to represent charged excitations of the insulator, we introduce bosonic operators pi and 
hi, representing the doubly-occupied and empty sites respectively. However, the singly- 
occupied site cannot be treated as a featureless vacuum, as we were able to in Section |T} 
Now we need a fermionic operator f a (the 'spinon') to represent the spin orientation of the 
singly-occupied site. For every site, we make the following correpondences for the four states 
in the Fock space 







tf|0> 


cm 






c\c\\0) 







It is now easy to verify that Eq. (2.7) is equivalent to the operator identification 

c a = (/i f +p)f a , 
provided we always project to states which obey the constraint 



(2.7) 



(2.8) 



(2.9) 



on every site. All physical observables are operators which stay within the subspace de- 



fined by Eq. (2.9): such operators are invariant under the following compact U(l) gauge 
transformation 



fa ^ fa 6 



p — > pe 



h -> he iC . 



(2.10) 



As we will show below, there will is an emergent gauge field _B„ in the effective theory 



9 



associated with this gauge transformation. The constraint in Eq. (2.9) will be the Gauss 



law of this gauge theory. These operator identities are related to those of the 'slave rotor' 
representation [8]. 

First, let us rewrite the Hubbard model in terms of these new bosonic and fermionic 



operators. The Hubbard Hamiltonian in Eq. (2.1) is now exactly equivalent to 



H[f,p,h] = -wJ^flfiaihi+plXht+pj) + H.c. 



£ (-Mrfft - h\K + 1) + U - [p\ Vl + h\K - ) , (2.11) 



provided our attention is restricted to the set of states which obey the constraint in Eq. (2.9) 



on every lattice site; note that the Hamiltonian in Eq. (2.11) commutes with constraints in 



(2.9), and so these can be consistently imposed. In the on-site terms Eq. (2.11) we have 



used the p and h operators to measure the electron density on each site 

We can now implement the commutation relations, the Hamiltonian, and the constraint 
in a coherent state path integral 



Z 



Vf ia (T)T>p i (T)Vh i (T)T>\ i (T) exp I - J drH[f,p, h] 



ftoTfiT + ^|r + + iX i(fLfia ~ PiPi + h\hi - 1) 



.dhj 



, (2.12) 



where the constraint in Eq. (2.9) is implemented using an auxilliary field Aj(r) which acts 
as a Lagrange multiplier. 



A key observation now is that the partition function in Eq. (2.12) is invariant under a 



site, i, and r-dependent U(l) gauge transformation 0( r ) where the fields transform as in 

(2.13) 



Eq. (2.10), and A transforms as 



A,- — y \j 



dr' 



In other words, A transforms like the temporal component of a U(l) gauge field. 

How do we obtain the spatial components of the gauge field? For this, we apply a 



"Hubbard-Stratonovich transformation" to the hopping term in Eq. (2.11). For this, we 



introduce another auxiliary complex field Qij{r) which lives on the links of the triangular 
lattice and replace the hopping term by 



E 

(ij) 



\Qa(r) 

w 



QuWfLfia - QvWhi+pbiht+pj) + H.c. 



(2.14) 



We now see from Eq. ( 2. 10 ) , that transforms under the gauge transformation in Eq. (2.10) 
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as 



Qij ^ Qij& 



(2.15) 



In other words, arg(Qij) is the needed spatial component of the compact U(l) gauge field. 
So far, we have apparently only succeeded in making our analysis of the Hubbard model in 



Eq. (2.1 ) more complicated. Instead of the functional integral of the single complex fermion 
Ci a , we now have a functional integral over the complex fermions fi a , the bosons Pi, hi, and 
the auxilliary fields Aj and Qij. How can this be helpful? The point, of course, is that the 
new variables help us access new phases and critical points which were inaccessible using the 
electron operators, and these phases have strong correlations which are far removed from 
those of weakly interacting electrons. 

The utility of the new representation is predicated on the assumption that the fluctuations 
in the auxiliary fields Qij and Aj are small along certain directions in parameter space. So 
let us proceed with this assumption, and describe the structure of the phases so obtained. 
We parameterize 



Qi 



-.Bi 



Xi = — iX — Bi 



(2.16) 



and ignore fluctuations in the complex numbers Q {j , and the real number A. With these 



definitions, it is clear from Eqs. (2.13) and (2.15) that Bij and B T form the spatial and 
temporal components of a U(l) gauge field, and so must enter into all physical quantities in 
a gauge invariant manner. The values of Q^ and A are determined by a suitable saddle-point 



analysis of the partition function, and ensure that the constraint (2.9) is obeyed. With these 



assumptions, the partition function separates into separate fermionic and bosonic degrees 
of freedom interacting via their coupling to a common U(l) gauge field (B iT ,Bij). In the 
continuum limit, the gauge fields become a conventional U(l) gauge field B^ = (B T ,B). 
The partition function of the gauge theory is 



Z 



VUa {r)V Pi {T)Vh t (r)VB iT (r)!?^ (r ) 



exp — / dr 



i 

U = XX (J^ - J - V + \ + %B ^) Pi + J2 h ^ +J + V + ^ - iB iT j hi 



Y,Q* j (h i +pl)e- iB v(h]+p j )+K.c. 

(ij) 



(2.17) 



Thus we have fermions fi a moving in a band structure which is roughly the same as that of 
the electrons in Eq. ( ]2.4[ ), the bosons pi, hi Hubbard-like Hamiltonian essentially identical in 
form to that in Section [l| and all particles are minimally coupled to a compact U(l) gauge 
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Vacuum \ 






Q/U 



FIG. 4: Possible phase diagram of the electron Hubbard model in Eq. (2.1) on the triangular 



lattice. This phase diagram is obtained by a mean-field analysis of the theory Cb in Eq. (2.17), 



similar to that for the boson Hubbard model in Chapter 9 of Ref. 4. We use the number operator 
n = p^p — h^h, which commutes with the boson hopping term, to characterize the Mott insulating 
stat es. O nly the Mott insulating lobes with n = —1,0,1 are compatible with the constraint in 
Eq. (2.9); these Mott insulating lobes have fermion density (fafa) = n + 1. 



field. 

We begin by neglecting the gauge fields, and computing the separate phase diagrams of 
Cf and Cb- 

The fermions are free, and so occupy the negative energy states determined by the chem- 
ical potential A. 

The phase diagram of Cb is more interesting: it involves strong interactions between the p 
and h bosons. It can be a analyzed in a manner similar to that of the boson Hubbard model 
(see Chapter 9 of Ref. HJ) , leading to the familiar "Mott lobe" structure shown in Fig. |1J 

At large values of Q/U we have the analog of the superfluid states of the boson Hubbard 



model, in which there is a condensate of the same operator as that in Eq. (1.9): 



?P=p + h j . (2.18) 

Note that ip is the ladder operator for the number operator n = p^p—h^h used to characterize 
the insulating phases in Fig. Q. The ip operator carries unit charge under the U(l) gauge 
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field (from Eq. (2.10), and so the superfluid phase, with (■0) 7^ 0, does not break any global 
symmetries (unlike the boson model of Section [T]). Instead it is a "Higgs" phase. In the 
presence of the Higgs condensate, the operator relation in Eq. (2.8 ) implies that c a ~ f a , and 
so the fermions carry the same quantum numbers as the physical electron. Consequently, 
the f a Fermi surface is simply an electron Fermi surface. Furthermore, the Higgs condensate 
quenches the fluctuations, and so there are no singular interactions between the Fermi 
surface excitations. This identifies the present phase as the familiar Fermi liquid, as noted 
in Fig. |4j We note that we can equally well identify this phase as a "confining" phase of the 
U(l) gauge theory, in which the ip boson has formed a bound state with the fermions, 
which is just the gauge-neutral c Q fermion. Indeed, as is well known, Higgs and confining 
phases are qualitatively the same when the Higgs condensate carries a fundamental gauge 
charge, as is the case here. 

Having reproduced a previously known phase of the Hubbard model in the U(l) gauge 
theory, let us now examine the new phases within the 'Mott lobes' of Fig. |4} In these states, 
the boson excitations are gapped, and number operator n = p'p—h'h has integer expectation 
values. The constraint in Eq. (2.9) implies that only n = —1, 0, 1 are acceptable values, and 
so only these values are shown. It is clear from the representation in Eq. (2.7) that any 
excitation involving change in electron number must involve a p or h excitation, and so the 
gap to the latter excitations implies a gap in excitations carrying non-zero electron number. 
This identifies the present phases as insulators. Thus the phase boundary out of the lobes 
in Fig. [4] is a metal-insulator transition. 

The three insulators in Fig. [4] have very different physical characteristics. 

Using the constraint in Eq. ( 2.9[ ) we see that the n = — 1 insulator has no f a fermions. 
Consequently this is just the trivial empty state of the Hubbard model, with no electrons. 

Similarly, we see that the n = 1 insulator has 2 f a fermions on each site. This is the just 
the fully-filled state of the Hubbard model, with all electronic states occupied. It is a band 
insulator. 

Finally, we turn to the most interesting insulator with n = 0. Now the electronic states 
are half-filled, with (flf a ) = 1- Thus there is an unpaired fermion on each site, and its 
spin is free to fluctuate. There is a non-trivial wavefunction in the spin sector, realizing an 
insulator which is a 'spin liquid'. In our present mean field theory, the spin wavefunction 
is specified by Fermi surface state of the fermions. Going beyond mean-field theory, we 
have to consider the fluctuations of the gauge field, and determine if they destabilize the 
spin liquid. The f a fermions carry the _B M gauge charge, and these fermions form a Fermi 
surface. The gapless fermionic excitations at the Fermi surface prevent the proliferation of 
monopoles in the compact U(l) gauge field: the low energy fermions suppress the tunneling 
event associated with global change in B^ fluxjHJ [TO]. Thus the emergent U(l) gauge field 
remains in a deconfined phase, and this spin liquid state is stable. These gapless gauge 
excitations have strong interactions with the f a fermions, and this leads to strong critical 
damping of the fermions at the Fermi surface which is described by a strongly-coupled field 
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theory [TTl - fT3] . The effect of the gauge fluctuations is also often expressed in terms of an 
improved trial wavefunction for the spin liquid [5]: we take the free fermion state of the f a 
fermions, and apply a projection operator which removes all components which violate the 



constraint in Eq. (2.9). This yields the 'Gutzwiller projected' state 



| spin liquid) = I ] J 



J] fl(k)fl(k) |0>, (2.19) 



\k<kp 



where the product over k is over all points inside the Fermi surface. 

Finally, we turn to an interesting quantum phase transition in Fig. |4} This is the transition 
between the spin liquid and the Fermi liquid at total electron density M = 1, which occurs 
at the tip of the n = Mott lobe. From the boson sector, this looks like a Higgs transition, 
of the condensation of a complex scalar ip as in Section [I], but in the presence of a fluctuating 
U(l) gauge field. However, the fermionic sector is crucial in determining the nature of this 
transition. Indeed, in the absence of the Fermi surface, this transition would not even exist 
beyond mean field theory: this is because the U(l) gauge field is compact, and the scalar 
carries unit charge, and so the confining and Higgs phases of this gauge theory are smoothly 
connected. So we have to combine the Higgs theory of a complex scalar with the gapless 
Fermi surface excitations. We can obtain the field theory for this metal-insulator transition 
by applying the methods of Section [T] to Eq. (2.17). The analog of the condition A p = Ah 
needed to obtain a density of one electron per site is now A + fi = 0. In this manner, we 
find the field theory [TJ [T3] 



C 



+ fi 



d_ 



e F 



+ *\ 



+ u\iP\ 4 + iB T M 



2m 



fa, 



(2.20) 



where the energy ep is to be adjusted to yield total fermion density M = 1. The transition 
is accessed by tuning s, and we move from a spin liquid with (ip) = for s > s c , to a 
Fermi liquid with (ip) ^ for s < s c . Within the the spin liquid phase, we can safely 



integrate out the gapped ip quanta, and so obtain the theory in Eq. (3.1) of the main 

s c have been studied [7J E], and an 



text. The critical properties of the theory at s 
interesting result is obtained: the Fermi surface excitations damp the gauge bosons so that 
they become ineffective in coupling to the critical b fluctuations. Consequently, the gauge 
bosons can be ignored in the ip fluctuations, and the transition is in the universality class of 
the 2+1 dimensional XY model. In other words, quite unexpectedly, the critical theory is the 
same as that of the superfluid-insulator transition of Section [Tj There are additional gapless 
excitations associated with the gauge field and the Fermi surface, but these are irrelevant 
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for the values of certain critical exponents. 
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